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Variational calculations for s-shell hypernuclei are performed by explicitly including E degrees 
of freedom. Four sets of YN interactions (SC97d(S), SC97e(S), SC97f(S) and SC89(S)) are used. 
The bound-state solution of ^He is obtained and a large energy expectation value of the tensor 
KN — T,N transition part is found. The internal energy of the *He subsystem is strongly affected 
by the presence of a A particle with the strong tensor AA'^ — EA'^ transition potential. 
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PACS numbers: 21.80.+a, 21.45.+V, 21.10.Dr, 13.75.Ev 

Few-body calculations for s-shell hypernuclei with 
mass number ^ = 3 — 5 are important not only to ex- 
plore exotic nuclear structure, including the strangeness 
degrees of freedom, but also to clarify the charac- 
teristic features of the hyperon-nucleon (YN) interac- 
tion. Although several interaction models have been 
proposed|jl|, ||, |j, the detailed properties (e.g. ^5*0 or 
^5*1 — ^ Di phase shift, strength of AN - SA^ couphng 
term) of the YN interaction are different among the mod- 
els. The observed separation energies (B/^'s) of light A 
hypernuclei are expected to provide important informa- 
tion on the YN interaction, because the relative strength 
of the spin-dependent term or of the AA^ — Y.N coupling 
term is affected from system to system. 

Recently, few-body studies for A — 3,4 hypernuclei 
have been conducted using modern YN interactions |4| 
^, |j. According to these developments, the Nijmegen 
soft core (NSC) model 97f (or 97e) seems to be compat- 
ible with the experimental B\s, though the calculated 



Ba for 



He* 



is actually slightly smaller than the 
experimental value. These few-body calculations, how- 
ever, have not yet reached a stage to calculate BaIaHg). 

If one constructs a phcnomenological central AA^ 
potential, which is consistent with the experimental 
Ba(3 H), BA(iH), BA(iHe), BaUB*) and SA(iHe*) val- 
ues as well as the Ap total cross section, that kind of 
potential would overestimate the BA{\tie) value [0, 1). 
This is known as an anomalously small binding of aHc. 
Though a suppression of the tensor forces |l^ or of 
the AA^ - SA^ coupling |l| was discussed to be a 
possible mechanism to resolve the anomaly, the prob- 
lem still remains an enigma JTsf due to the difficulty 
of performing a complete five-body treatment. Only 
one attempt was made, using a variational Monte Carlo 
calculationQ with the NSC89 YN interaction. Though 
NSC89 well reproduces both the experimental i?A(AH)@] 
and i3A(AH)[§[ ^ values as well as the experimental Ap 
total cross section, a bound-state solution of ^He was not 
found. In view of the aim to pin down a reliable YN in- 
teraction, a systematic study for all s-shell hypernuclei 



is desirable. 

The A'^A^ tensor interaction due to a one-pion-exchange 
mechanism is the most important ingredient for the bind- 
ing mechanisms of light nuclei. More than a third, or 
about one half, of the interaction energy comes from the 
tensor force for the "^HeljlSl, Since the pion-(or 

kaon-) exchange also induces the AA^ — J^N transition 
for the YN sector, both the A^A^ and AA^ - T,N tensor 
interactions may also play important roles for light hy- 
pernuclei. If this is the case, the structure of the core 
nucleus (e.g. '^He) in the hypernucleus (aHc) would be 
strongly influenced by the presence of a A particle. 

The purpose of this letter is twofold: First is to perform 
an ab initio calculation for ^He as well as A = 3, 4 hyper- 
nuclei explicitly including S degrees of freedom. Second 
is to discuss the structural aspects of ^He with an appro- 
priate YN interaction which is consistent with all of the 
s-shell hypernuclear data. 

The Hamiltonian (H) of a system comprising nucleons 
and a hyperon (A or E) is given by 2 x 2 components as 



H 



Ha ^s-a 
Va-s 



(1) 



where Ha{H^) operates on the A-(I]-)component and 



Va-j 



1=1 



(2) 



We employ the G3RS potential|]T8| for the A^A^ interac- 
tion and the SC97d(S), SC97e(S), SC97f(S) or SC89(S) 
potential[|9) for the YN interaction, where all interac- 
tions have tensor and spin-orbit components in addition 
to the central one. We omit small nonstatic correction 
terms {{L ■ S f and terms) in the G3RS A^A^^ interac- 
tion and odd partial- wave components in each interaction 
in order to focus on the main part of the interaction in 
the even parity state. The calculated binding energies 
for hght nuclei (^H, ^H, ^He and ^He) are 2.28, 7.63, 
6.98 and 24.57 MeV, respectively. The YN interactions 
have Gaussian form factors whose parameters are set to 
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reproduce the low-energy S matrix of the corresponding 
original Nijmegen YN interactions [po]. These Gaussian 
form factors help to save significant computer time. 

The binding energies of various systems are calculated 
in a complete A-hody treatment. The variational trial 
function must be flexible enough to incorporate both the 
explicit E degrees of freedom and higher orbital angular 
momenta. The trial function is given by a combination 
of basis functions: 

*.7j\/ta/t = Z]^=i Cfc</?fc, with 

ipk=A{G{x; Ak)[9Lk{x; Uk, Kk)xsAjM'nkTMT] ■ 

Here, A is an antisymmetrizer acting on nucleons and 
XSfc {VkTMr) is the spin (isospin) function. r]kTMT has 
two components: upper (lower) component refers to the 
A-(E-)component. The abbreviation x = {xi, • • ■ , xa-i) 
is a set of relative coordinates. A set of linear variational 
parameters (ci, • • ■ , cn) is determined by the Ritz varia- 
tional principle. 

A spatial part of the basis function is constructed by 
the correlated Gaussian(CG) multiplied by the orbital 
angular momentum part Ol{x), expressed by the global 
vector representation(GVR) pi| . CG is defined by 

= exp I - 2 Y.„3=i {Ak),j Xi ■ Xjj. 

The (A — l) X (A — l) symmetric matrix (Ak) is uniquely 
determined in terms of the intcrparticle correlation pa- 
rameter (akij)- The GVR of 9Lk{x;uk, Kk) takes the 
form 

9L,{x,Uk,Kk)^vl''''+''''YL,{vk), with 

The Ak and Uk are sets of nonlinear parameters which 
characterize the spatial part of the basis function. Al- 
lowing the factor v'^^'' {Kk ^ 0) is useful to improve the 
short-range behavior of the trial function. The value of 
Kk is assumed to take or 1. The variational parame- 
ters are optimized by a stochastic procedure. The above 
form of the trial function gives accurate solutions. The 
reader is referred to Rcf. jl^, ^ for details and recent 
applications. For the spin and isospin parts, all possible 
configurations are taken into account. 

Table | lists the results of the A separation energies. 
The scattering lengths of the ^So{as) and ^Si{at) states 
for each YN interaction are also listed in Table |[ where 
the interactions are given in increasing order of \as\ (and 
in decreasing order of \at\)- The SC89(S) interaction pro- 
duces no or very weakly bound state for , aHc* or 
^He. For the SC97d-f(S), the BA(^He) value is about 
2 — 3 MeV. This is a first ab initio calculation to pro- 
duce the bound state of ^He with explicit E degrees of 
freedom. 
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FIG. 1: Density distributions of A'^, A and E for |^He as a 
function of r, the distance from the center-of-mass of *He. 
The SC97e(S) YN interaction is used. Note that the E- 
distribution has been multiplied by a factor of 10 to clarify 
the behavior. The dotted line is taken from Ref. [|22|. 



The order of the spin doublet structure of the A — 4 
system is correctly reproduced for all YN interactions; 
the ground (excited) state has spin-parity, J"^ = 0+(l"'") 
for both isodoublet hypernuclei and Xlie. Although 
the strengths of the ^5*0 and ^Si interactions of the 
SC97d(S) are almost the same as each other, the energy- 
level of the 0+ state is clearly lower than that of the 1+ 
state. All of the ^4 = 3 bound states given in Table || have 
J'^ = in agreement with experiment. No other bound 
state has been obtained for all of the YN interactions. 
For the SC97e(S), the differences between the calculated 
and experimental B\ values are the smallest among the 
YN interactions employed in the present study. 

Table || lists the probability, (in percentage), of 
finding a E particle in the system. The sizable amount 
of Ps(AHe)'s is obtained. This implies that the A — E 
coupling plays an important role, even for the de- 
spite a large excitation energy of the core nucleus, *He 
(with the isospin 1), in the E-component. For the A = 4 
system, the Ps's of the 0+ state are about 1 — 2%, except 
for the SC89(S), while the Ps's of the 1+ state are nearly 
equal to or smaller than that of the 0+ state. 

Figure |^ displays the density distributions for ^He us- 
ing SC97e(S), and of N, A and E from the center-of-mass 
(CM) of '^He. Figure |^ also shows the A-distribution 
obtained from the Isle A — a potential |2^. The ex- 
perimental pionic decay width of ^He suggests that the 
A-distribution should spread over a rather outer region 
compared to the distribution of the a, as was discussed 
in Ref. [Q. The present curve of the A-distribution is 
similar to that obtained by the Isle potential. The E- 
distribution has a shape similar to the A^-distribution. 
The root-mean-square (rms) radii of N^ A and E from 
the CM of the "^He are 1.5, 2.9 and 1.6 fm, respectively. 
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TABLE I: A separation energies, given in units of MeV, of yl = 3 — 5 A-hypernuclei for different YN interactions. The scattering 
lengths, given in units of fm, of ^5'o(aa) and 5*1 (at) states are also listed. 



YN 


as 


at 


SA(iH) 


BA(iH) 




BA(iHe) 


SA(iHe*) 


BA(iHe) 


SC97d(S) 


-1.92 


-1.96 


0.01 


1.67 


1.20 


1.62 


1.17 


3.17 


SC97e(S) 


-2.37 


-1.83 


0.10 


2.06 


0.92 


2.02 


0.90 


2.75 


SC97f(S) 


-2.82 


-1.72 


0.18 


2.16 


0.63 


2.11 


0.62 


2.10 


SC89(S) 


-3.39 


-1.38 


0.37 


2.55 


unbound 


2.47 


unbound 


0.35 


Expt. 






0.13 ±0.05 


2.04 ±0.04 


1.00 ±0.04 


2.39 ±0.03 


1.24 ±0.04 


3.12 ±0.02 



TABLE IL Probabilities, given in percentage, of finding a E particle in A = 3 — 5 A-hypernuclei for different YN interactions. 



YN 








ft(iHe) 




ftdHe) 


SC97d(S) 


0.06 


1.27 


1.37 


1.24 


1.35 


2.04 


SC97e(S) 


0.15 


1.49 


0.98 


1.45 


0.96 


1.55 


SC97f(S) 


0.23 


1.88 


1.09 


1.83 


1.08 


1.87 


SC89(S) 


0.65 


3.73 


unbound 


3.59 


unbound 


1.33 



Table III lists the energy expectation values of the ki- 
netic and potential energy terms for A-He. The contribu- 
tions from the spin-orbit and the Coulomb potentials are 
not shown in the table, though the calculations include 
them. Here, Tc is the kinetic energy of the core nucleus 
(c) subtracted by the CM energy of c: 



Tr. = 



A-l 

E 

1=1 



A-l 
= 1 Pi 



2mis! 2{A — Vjmisi 



(6) 



The kinetic energy of the relative motion between the Y 
and the CM of c is given by 



Tv-r. — 



2^Y 



+ {'mY —mtCjc^ 



(7) 



{A—l)raj^raY 



is the reduced mass for the 



where iiy — rn — ri ; 

Y+c system and tvy-c is the canonical momentum of the 
relative coordinate between Y and c {Y — A, S). Ty-c 
also counts the difference in the rest-mass energy between 
A and E. Each potential part {V) takes account of a 
summation over appropriate particle pairs (see Eq. (|^) 
for example). The energy expectation values of the first 
three columns in Table III are written as 



(^'a|0|*a)±(*e|0|*s> 



(8) 



where the upper (lower) component of the ^ jmtMt 
denoted by \['a ('I'e)- The first (second) term of each 
element ((Tc), (Ty-c) or (Vnn)) in Table III represents 
the first (second) term of Eq. (^). The energy of the ''He 
subsystem changes a lot from that of the isolated one. 



AEr 



[{Tc) ± (Vnn] 
4.7MeV. 



iHe 



{Tc) + {Vnn)), 



He 



(9) 



This difference is considerably large despite the fact that 
the rms radius of N from the CM of '^He for the \lle 
hardly changes from that for "'He. (Both radii are 1.5 fm.) 
Most of the change is due to a reduction of the energy 
expectation value of the tensor NN interaction. 



((l^ArAr(tensor))^5^^ - (^(VwA' (tensor)} 



'He 



2.9MeV. 



On the other hand, the tensor AA^ — J^N transition part 
has a surprisingly large energy expectation value (about 
—20 MeV). This large couphng energy makes ^He bound 
in spite of both the energy loss of AEc and the extremely 
high energy of the E-component ( ^^^^ 600 MeV). 

The calculated wave function is divided into orthog- 
onal components according to the total orbital angular 
momentum (i), the total spin (S*), the core nucleus spin 
(Sc) and the core nucleus isospin (Tc). Table displays 
the probability of each component for \ile. The table 
also lists the probability of S'-state or of D-state for '*He. 
The sizable amount of probability of the E-component is 
found in the £)-state while the sum of IJ-state probabili- 
ties in the A-component is slightly smaller than that for 
"'He. Moreover, though the presence of a A in ''He with 
the strong tensor AN—'SN transition potential influences 
the structure of the -D-state component and reduces the 
energy expectation value of the tensor NN interaction, 
the large coupling energy (Va-e) of the tensor part bears 
the bound state of ^He instead. 

In summary, we have made a systematic study of all 
s-shell hypernuclei based on ab initio calculations us- 
ing YN interactions with an explicit S admixture. The 
bound-state solution of ^He was obtained. As the Ref. ||] 
claimed, though there is none of the interaction models 
to describe very precisely the experimental i?A's, the five- 
body calculation convinced us that the anomalous bind- 
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TABLE III: Energy expectation values of the kinetic and potential energy terms for aHc, given in units of MeV. The SC97e(S) 
YN interaction is used. For each potential part, a summation over appropriate particle pairs is taken into account (see Eq. ^ 
for example) and two central (^E and '^E) and a tensor (^E) components are listed separately. The first (second) term of each 
element {Tc), (Ty-c) or {Vnn) represents the first (second) term of Eq. (O = Tc,Ty-c or Vivjv). The energy expectation 
values of (Tc) and three (Vivjv)'s for isolated *He are 84.86, -33.22, -33.05 and -43.93 MeV, respectively. 



(Tc) (Ty-c) 


(Vjviv) 


{Vna) 


2{Va-s) 


{Vni:) 






83.43+2.74 9.11 + 3.88 


-33.14-0.35 


-3.97 


-0.02 


0.07 


(Central, 


'E) 




-32.03 - 0.27 


2.98 


-1.02 


1.56 


(Central, 






-40.91 - 0.12 


-2.24 


-19.51 


0.87 


(Tensor, 


'E) 



TABLE IV: Probabihty, given in percentage, of each compo- 
nent with the total orbital angular momentum (L) , total spin 
(5), core nucleus spin (Sc) and core nucleus isospin (Tc) in A- 
or in E-component for \VLe. The SC97e(S) YN interaction 
is used. The probability in S- or in _D-state for *IIe is also 
listed. 









L 


= 




L = 2 








S 


_ 1 

^ 2 


s 


3 Q 5 








=0 


Sc = 1 


Sc = 1 


Sc 2 Sc 2 


iHe 














{T.= 


0)(} 


5A 


89.14 


0.03 


0.19 


3.74 5.36 


(T,= 






0.10 


0.09 


1.34 


~ 0.01 


^He 






89.56 






10.44 



ing problem would be resolved by taking account of the 
explicit S admixture. The Bj^ values, obtained by us- 
ing SC97e(S), are the closest to the experimental values, 
among the YN interactions employed in this study. A 
sizable amount of Pe was obtained, even for the ^He, in 
spite of the large excitation energy of '^He. The contribu- 
tion of the energy from the tensor A7V — YiN coupling is 
quite large, and this coupling is considerably important 
to make ^He bound. This is novel finding in contrast 
with the Brueckner-Hartree-Fock calculation [|l9j. The 
present study for ^He is a first step toward a detailed 
description of light strange nuclear systems. The core 
nucleus, ^He, is no longer rigid in interacting with a A 
particle. A similar situation can occur for strangeness 
S — —2 systems. Investigations into the strength of the 
AA interaction based on the experimental data of the 

6t 



binding energy for double A hypernuclei (e.g. aaHb |23 ) 
should take account of the energy reduction of the core 
nucleus {/S.Ec for ^AHe is expected to be larger than the 
present Ai^c for ^He). 
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